Andreev Peaks and Massive Magnons in Cuprate SNS junctions 
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The projected SO(5) theory (pSO(5)) is used to resolve the puzzle of two distinct energy gaps 
in high T c Superconductor-Normal-Superconductor junctions. Counter to conventional theory of 
multiple Andreev reflections (MAR), the differential resistance peaks are associated with the anti- 
ferromagnetic resonance observed in neutron scattering, and not with Cooper pair breaking. The 
pSO(5) and MAR theories differ by the expected tunneling charges at the peaks. We propose that 
shot noise experiments could discriminate against the conventional interpretation. 
PACS numbers: 74.20.-z,74.65.+n 
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In current transport through high T c superconductor 
junctions, there seem to be two energy scalesa. The upper 
energy is seen in tunneling conductancecl, and is identi- 
fied with the "pseudogap" A p which appears in magnetic 
resonances and photoemmissiona. A lower gap, which 
scales differently with hole doping, manifests as peaks 
in the differential resistance of Superconductirur-Normal- 
Superconducting (SNS) Josephson junctional. These 
peaks have been interpreted using the conventional the- 
ory of multiple Andreev reflections, following Klapwijk, 
Blonder, and Tinkham (KBT)B. 

KBT theory treats two conventional superconductors 
with a single s-wave BCS quasiparticle gap A, separated 
by a free electron metal. Electrons traversing the metal 
are Andreev reflected back as holes, gaining energy incre- 
ments eV at each traversal (as depicted in Fig. |l|). Peaks 
in the differential resistance appear at voltages 2A/ne, 
and are due to the (E — A) -1 / 2 singularity in the quasi- 
particles' density of states. 
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FIG. 1. KBT Theory. Differential resistance peaks of 
n = 6 (left diagram), and n = 5 (right diagram), involve a 
cascade of n Andreev reflected charges traversing the normal 
metal. Singular dissipation is due to emission of quasiparticles 
above the s-wave gap. Filled (empty) circles denote electrons 
(holes) in the normal barrier. 



However, in cuprate SNS junctions, such as YBa 2 Cu3- 
6 .6 - YBa2Cu2.5sfeo.45Oy -YBa 2 Cu 3 06.6 examined by 
Nesher and KorenQ , application of KBT theory is prob- 
lematic. A naive fit to KBT expression faces the two gaps 



puzzle, i.e. an "Andreev gap" is of order A « 16mcV, 
while the tunneling gap is about three times largerQ, and 
scales differently with T c . Without perfect alignment 
of the interfaces, it, is hard to understand the observed 
sharpness of peaksu since the d-wave gap is modulated 
for different directions. Moreover, the barrier is by no 
means a "normal" metal devoid of interactions: it is an 
underdoped cuprate with antiferromagnetic correlations 
and strong pairing interactions as evidenced by a large 
proximity effecttl. 

The purpose of this Letter is to provide an alternative 
explanation for the differential resistance peaks seriesO, 
which takes into account the strong correlations in the 
pseudogap regime. Our analysis resolves the two energy 
scales puzzle. ^ 
We employ the projected SO(5) (pSO(5)) modeel, which 
is a strong coupling effective Hamiltonian. It describes 
the dynamics and interactions of four primary bosonic 
modes of cuprates: preformed hole pairs and massive spin- 
one magnons. 

A differential resistance peaks series is found at bias volt- 
ages V n = A s /(en), n = 1,2,... where A s is the an- 
tiferromagnetic resonance energy. This resonance has 
been directly measured by inelastic neutron scattering. 
The peaks are thus associated with emmission of magnon 
pairs at the resonance threshold, and not with pair break- 
ing, as in KBT theory. We note that other predictions to 
observe magnons (also called 7r-modes) in various cuprate 
junctions were madeE3iij, but await experimental con- 
firmation. We propose that measurement of the excess 
shot noise below the peaks, could discriminate against 
the latter interpretation. pSO(5) theory predicts tunnel- 
ing charge 2ne below the n-th peak, while KBT theory 
expects charge ne. 

Degrees of freedom: At energies below the pseudogap A p , 
preformed hole pairs (with internal d-wave symmetry), 
describe the primary charge degrees of freedom in the un- 
derdoped regimeta. The hole pairs are bosons, and their 
phase fluctuations are controlled by the two dimensional 
superconducting stiffness p c , as measured by the Lon- 
don penetration depth. At T c , the pairs Bose condense 
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and long range phase coherence is established. This sce- 
nario can explained the empiricai*elations T c oc p c , which 
have been observed in cupratesli3 at low doping concen- 
trations. The other low energy charge excitations are 
fcrmionic quasiparticles near the d-wave nodes. These 
have a smooth density of states which decreases below 
A p . 

Additional bosonic excitations below the pseudogap en- 
ergy scale, are antiferromagnetic spin fluctuations i.e. 
magnons. Massive spin one magnons have been observed 
in inelastic neutron scattering^ in YBa2Cu30g + 5. They 
manifest as a sharp resonance in the spin correlation func- 
tion S aa >, which near the antiferromagnetic wavevector 
q w 7? has the form 

uj z — cr(q — -k) 1 — A z s 

Here c is the spin wave velocity, and s is a normal- 
ization factor. iT-he doping dependent resonance energy 
A s (5) increases^ between A s (0.5) = 25meVp-(with T c = 
52°K), and A s (l) = 40meV, (at T c = 92° Itfi. The pro- 
jected SO (5) theory. The large onsite Hubbard repulsion 
between electrons is imposed by an apriori projection of 
doubly occupied states from the Hilbert spacec3. 
The undoped vacuum |0) is a half filled Mott insulator 
in a quantum spin liquid state. The pSO(5) vacuum 
possesses short range antiferromagnetic correlations. A 
translationally invariant realization of |0) on the micro- 
scopic square.lattice, is the short range resonating valence 
bonds stateBEl 

Out of this undoped vacuum, b h create charge 2e bosons 
(hole pairs) with internal d-wave symmetry under rota- 
tions, and }r a , a — x, y, z create a triplet of antiferromag- 
netic, spin one magnons. 
The lattice pSO(5) Hamiltonian is 

ncharge = (£c _ 2/l) ^ biAi _j*y j (b) llbhl + h. c .) 

« (ij) 
ia a(ij) 



n mt = wY / --[b{ l b hl + ^2blb 




(2) 

where : () : denotes normal ordering, and nf — (b ia + 
b ia )/ ' \[2 is the Neel spin field. TL mt describes short range 
interactions between bosons, and 7i ul describes the 
long range Coulomb interactions. H^ erm describes cou- 
pling to the nodal (fermionic) quasiparticles, which con- 
tribute to a large, but smooth, conductance background. 
Here we will concentrate on the conductance singulari- 
ties, and will not compute the fermionic background. 



The mean field approximation to Eq. (Q) is 
straightforward!!. It amounts to replacing b*- — > (6^), 
7 = h,a, and minimizing U charge + H spm + H mt with 
respect to (&*,,). There is a first order transition between 
the two primary mean field phases on the square lattice 
at p = p c , where 

ptc = ^(e c - e s ) - (J c - 2J S ), (3) 

p c is of the order of the Hubbard interaction scale. 
At p < p c we have an undoped Mott insulator with no 
hole pair bosons, and where the magnons Bose-condense. 
The condensate supports a finite staggered magnetiza- 
tion 

\(n a )\ 2 = (2J s -±e s )/W p<p c (4) 

There are two linear spin wave modes u = c|q|, with 
where c = \^2J s /h is the semiclassical spinwave velocity 
of the Heisenberg antiferromagnet. 

At p > p c the ground state becomes doped with hole 
pairs which Bose-condense into a superconducting phase 
with an order parameter 

| (b\) | 2 = ( Ji - p c + 2 J s - e s /2)/W p> p c (5) 

Long range interactions in T-i Coul j frustrate the first order 
transition ancLcreate intermediate (possibly incommen- 
surate) phaseso, which we shall not discuss here. 
The mean field phase stiffness is given by p c — J c (b\} 2 , 
and therefore Eq. (|^) explains why p c increases with 
chemical potential (and doping) in the underdopsd su- 
perconducting regime, as observed experimentally^. 
Analysis of the linear quantum fluctuations about mean 
field theory determines the magnon dispersion i.e. the 
poles of Eq. (|l|) . The mean field magnon gap is found to 
be 

A s = 2y/(fi- At c ) (/i- Pc + ^Js) (6) 

which by Eq. (g) implies that A 2 oc p Cl T c . Thus the 
pSO(5) mean field theory can explain the systematic in- 
crease of A s with T c which is observed by Fong et. al&3. 

The cuprate SNS junction. We consider a junction, where 
the barrier (N) has no superconducting or magnetic or- 
der (b{) = 0, (n a ) = 0. We derive on general grounds the 
form of the effective tunneling Hamiltonian between su- 
perconductors as follows. An integration of the barrier's 
charged bosons bh out of the path integral results in an ef- 
fective action S tun which couples the charges of the two 
superconductors . S tun [bh L ,bh R ,b a ] explicitly depends 
on the hole pairs bosons on the left and right interfaces, 
and on the magnons in the barrier. By charge conserva- 
tion, an expansion of S tun as a power series leaves only 
terms with equal number of bhS and b h 's. By spin con- 
servation, the magnon terms are singlets and hence at 
least bilinear in n a . 
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FIG. 2. pSO(5) theory for Andreev peaks in cuprate 
SNS junctions. Three hole pairs co-tunneling from left to 
right, generate a pair of magnons. At the antiferromagnetic 
resonance threshold 6eV = 2A S , this process contributes to 
the n = 3 peak of the differential resistance. The diagram 
contains lowest order contributions of hole pairs-magnon in- 
teractions to the tunneling vertex T3. 



This expansion leads to a series of tunneling terms. For 
the Andreev peaks we retain only the leading order terms 
(in &t t fr) which are 



U tun-mag = _Y^(A n + Ji) 
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Jv n describes a simultaneous tunneling of n hole pairs 
from the left to the right superconductor, coupled to a 
magnon pair excitation. T n is the tunneling vertex func- 
tion, which depends on the bosons positions. 
The energy transfer mechanism is depicted diagrammat- 
ically in Fig.^. We do not compute T„'s which depend on 
the details of the barrier and the interfaces. A "good" N 
barrier is defined to have sizeable T n , if multiple pair tun- 
neling terms are to be observed. This requires a thin bs 
rier with slowly decaying spin and charge correlations 
It is important to note that multiple pair tunneling, i.e. 
the differential resistance peaks at n > 1, depends on 
strong anharmonic interactions between the hole pairs 
and magnons. These interactions are an essential part 
of the pSO(5) theory as modelled by H int in Eq. ^. 
The juBGjtion's conductance is calculated in the standard 
fashionlia: the bias voltage V transforms the left bosons 
bh L — > e z2eVt bh L , which yields time dependent operators 
An(t). The current is calculated by second order pertur- 
bation theory in J{ tun ~ ma 9 yielding 



I = J22neX r n et (2eV) 



rH=* dte^{[A n (t),Al}) (8) 
Jo 



For singular contributions / sm 9 ; we ignore supercon- 
ducting condensate fluctuations b* h — (6^), which have 
a smooth spectrum. Similarly, we ignore the frequency 
dependence of T n (uj). Setting b* R — > (6^) 
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S(q,iu + 2neV + iO + )S(-q,iu) (9) 



where the barrier dimensions are d x W (see Fig.|), and 
^ w is a Matsubara sum. 

For a nearly antiferromagnetic "N" barrier, T„(x — x') 
in (0) decays slowly with the distance between magnons. 
Thus for a narrow barrier d « W, the magnons are 
excited at q y ss 0, and the momentum sum reduces to a 
one dimensional sum over q x . At zero temperature we 
obtain 

F in B = Y^2ne{bl) in \T n [Q]\ 2 

n 

dq x S(2neV - 2y/c 2 ql + Af) 
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The last expression emphasizes the singular form of 
I smg (V,A s ) at the peaks. For a large background con- 
ductance dl/dV » dl smg /dV, the inverse square root 
singularities in I sm 3 create peaks in the differential re- 
sistance dV/dl at voltages 



V n = A s /(ne), n=l,2, 



Q n = 2ne (11) 



where Q n is the excess tunneling charge below the n-th 
peak. Note that Q n changes in increments of 2e. The 
differential resistance peak series is depicted in Fig. ||, 
for weak broadening of the singularities and an arbitrary 
set of coefficients t n . 
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n < 5, and a background con- 



ductance of unity. Below the n-th peak, the excess tunneling 
charge is 2ne, rather than BTK's ne. 
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Discussion. We have seen that magnon pair creation in- 
duces peaks in the differential resistance which are similar 
in appearance to the Andreev peaks of the KBT mech- 
anism. The crucial difference is that here the singular 
dissipative process does not involve Cooper pair breaking, 
but low energy antiferromagnetic excitations. 
In KBT theory for two identical superconductors, the 
peaks appear at voltages V^ 37 = 2A/(ne), n = 
1,2,... which are the upper threshold for tunneling of 
charges Q n = ne. Thus, KBT allows both even and odd 
number of electron charges to participate in the multiple 
Andreev reflection process, as depicted in Fig. [|. These 
charges change in increments of e at each peak. There- 
fore a decisive discrimination between the processes of 
Fig.jl] and Fig.|] would be measurements of the excess 
tunneling charge increments at the peaks. A feasible 
method would perhaps be low temperature shot noiseH S, 
which measures the tunneling charges via the relationE2l 
S = 2Q n I(V n ). We eagerly look forward to the results of 
such experiments. 
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